1 



Parton Distribution Functions 

Davison E. Soper a * 

institute of Theoretical Science, University of Oregon, 
Eugene, Oregon 97403, USA 



Talk at Lattice 96 Conference, St. Louis, June 1996. 

Parton distribution functions give the probability to find partons (quarks and gluons) in a hadron as a function 
of the fraction x of the proton's momentum carried by the parton. They are conventionally defined in terms of 
matrix elements of certain operators. They are determined from experimental results on short distance scattering 
of the partons. Integrals of these functions weighted with x" are calculable using lattice QCD. Some simple 
models for their behavior have implications that could also be tested in lattice QCD. 



1. Intuitive meaning of parton distribution 
functions 

Let da be a cross section involving short dis- 
tances. For instance, we may consider the process 
hadron A + hadron B — > jet + X at the Fermilab 
collider. Let the jet have a high transverse mo- 
mentum Pt- Intuitively, the observed jet begins 
as a single quark or gluon that emerges from a 
parton-parton scattering event with large Pt, as 
illustrated in Fig. [p. (Typically, this parton re- 
coils against a single parton that carries the op- 
posite Pt-) The large Pt parton fragments into 
the observed jet of hadrons. 

The physical picture illustrated in Fig. [l] sug- 
gests how we may write the cross section to pro- 
duce the jet as a product of three factors. A par- 
ton of type a comes from from a hadron of type 
A. It carries a fraction xa of the hadron's mo- 
mentum. The probability to find it is given by 
f a /A( x A) dxA- A second parton of type b comes 
from a hadron of type B. It carries a fraction xb 
of the hadron's momentum. The probability to 
find it is fb/ b{xb) dxB- The functions f a /A{%) 
are the parton distribution functions that are the 
subject of this talk. The third factor is the cross 
section for the partons to make the observed jet, 
da. This parton level cross section is calculated 
using perturbative QCD. 
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Figure 1. Hadron A + hadron B — > 2 partons. 



1.1. Factorization 

We have been led by the intuitive parton pic- 
ture of Fig. |l| to write the cross section for jet 
production in the following form 



da 

dPr~ ~ 

^2 Jdx A f a /A(xA,n) Jdx B fb/B(%B,H) 



(i) 



da 
dPr~' 



Here the parton level cross section has a well be- 
haved expansion in powers of a s , 



da 
dPr 



A 



N 



^ ( ~~~ 1 H n (xa,x b ,Pt;cl, b; fi). (2) 
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Figure 2. Factorization for hadron collisions. 



The coefficients Hn are calculable in perturbative 
QCD. 

The principle of factorization asserts that 
Eq. (0) holds up to corrections of order 

• (m/Pr) n where to is a typical hadronic mass 
scale and the power n depends on the 
process, and 

• (a s (/j,)) L from truncating the expansion of 
da/dP T - 

For our purposes, we can regard factorization as 
an established theorem of QCD, although this 
subject is not without its loose ends. A review 
may be found in Ref. 0. 

As we have seen, Eq. ([!]) has a simple intuitive 
meaning. However, the appearance of a param- 
eter /i in Eq. (Q) hints that there is more to the 
equation than just a model. The parameter [i, 
which has dimensions of mass, is related to the 
renormalization of the strong coupling a s (/i) and 
of the operators in the definition of the parton dis- 
tribution functions f a /A( x A, A*)- (Often, one uses 
two separate parameters in these two places.) 

At the Born level, the parton level cross section 
da 1 1 cLPt is calculated in a straightforward man- 
ner. At the next-to-leading order and beyond, 
the calculation is not so straightforward. Various 
divergences appear in a naive calculation. The di- 
vergences are removed and the dependence on the 
scale [i appears in their place. The precise rules 
for calculating da/dPr follow once you have set 
the definition of the parton distribution functions 
fa/Ai x i f)- These rules enable one to do practical 



calculations. For example, for jet production, the 
first two terms in da/dPx are known in the form 
of computer code that puts together the pieces 
of Eq. (|l|) and produces a numerical cross section 
§• 

1.2. Reality Check 

Sets of parton distribution functions, one func- 
tion for each kind of parton in a proton, are pro- 
duced to fit experiments. We will examine the 
fitting process in a later section. Fig. || is a graph 
of the gluon distribution and the up-quark distri- 
bution in a proton, according to a parton distri- 
bution set designated CTEQ3M [|. The figure f§ 
shows x 2 f a / A (x, /j,) for a — g and a — u, A = p. 
Note that 

J dx xf a/A {x,n) = J rflogx x 2 f a / A (x,n), (3) 

so the area under the curve is the momentum frac- 
tion carried by partons of species a. 
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Figure 3. Gluon and up quark distributions in the 
proton according to the CTEQ3M parton distri- 
bution set. 



1.3. Significance 

Knowledge of parton distribution functions is 
necessary for the description of hard processes 
with one or two hadrons in the initial state. With 



3 



two hadrons in the initial state, as at Fermilab or 
the future Large Hadron Collider, observed short 
distance cross sections take the form 

da ~ (4) 
^2 J dx A fa/A(xA, /-0 J dx B h/B(xB,n) da. 

a.b 

With one hadron in initial state, as in deeply in- 
elastic lepton scattering at HERA (Fig. ||), the 
cross section has the form 

da~^2 dxA f a /A(xA, /x) da. (5) 

In either case, one has no predictions without 
knowledge of the parton distribution functions. 




Figure 4. Deeply inelastic scattering. 

The essence of Eqs. ([!]) and (^) above is that in 
high energy, short distance collisions a hard scat- 
tering probes the system quickly, while the strong 
binding forces act slowly. Thus one needs to 
know probabilities to find partons in a fast mov- 
ing hadron as seen by an approximately instan- 
taneous probe. This is the information encoded 
in the parton distribution functions. It should be 
evident that this information is not only useful, 
but also valuable in understanding hadron struc- 
ture. Parton distribution functions are not every- 
thing, however. They provide a relativistic view 
only, a view quite different from the view that 
might be most economical for the description of 
a hadron at rest. Furthermore, they provide no 
information on correlations among the partons. 



2. Translation to operators 

We are now ready to examine the technical defi- 
nition of parton distribution functions. There are, 
in fact, two definitions in current use. I will de- 
scribe the MS definition, which is the most com- 
monly used. There is also a DIS definition, in 
which deeply inelastic scattering plays a privi- 
leged role. The interested reader may consult the 
CTEQ Collaboration's Handbook of Perturbative 
QCD H for information on the DIS definition. 
There are also different ways to think about the 
MS definition. In this talk, I define the parton 
distributions directly in terms of field operators 
along a light-like line, as in An equivalent 
construction from a different point of view may 
be found in 0. As we will see, moments of the 
parton distribution functions are related to ma- 
trix elements of certain local operators, which ap- 
pear in the operator product expansion for deeply 
inelastic scattering. This relation could be also 
used as the definition. 

I distinguish between the parton distribution 
functions f a /Ai x i^) & n d the structure functions 
Fi(a;,Q 2 ), F 2 ( Xj Q 2 ), and F 3 (x,Q 2 ) that are mea- 
sured in deeply inelastic lepton scattering. 

2.1. Null coordinates 

I will use null plane coordinates and momenta 
defined by 

x ± = (x°±x 3 )/V2 1 P ± = (P°±P 3 )/V2. (6) 

Imagine a proton with a big P + , a small P~ , 
and Pt = 0. The partons in such a proton move 
roughly parallel to the x + axis, as illustrated in 
Fig. [|. One can treat x + as "time," so that the 
system propagates from one plane of equal x + to 
another. For our fast moving proton, the interval 
in x + between successive interactions among the 
partons is typically large. 

Notice that the invariant dot product between 
P^ and x^ is 

P ■ x = P + x~ +P~x + - P T ■ x T ■ (7) 
Thus the generator of "time" translations is P~ . 

2.2. Null plane field theory 

We will define parton distribution functions by 
taking a snapshot of the proton on a plane of 
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Figure 5. World lines of the partons in a fast- 
moving proton. 



equal x + in Fig. |^. To motivate the definition, 
we use field theory quantized on planes of equal 
x + ||. This quantization uses the gauge A + = 0. 
Then the unrenormalized quark field operator tpo 
is expanded in terms of 

• quark destruction operators b and 

• antiquark creation operators d) 

using simple spinors w(s) normalized to w^w = 1: 



±7 1 + ip {x + ,x ,x T ) 

i r 00 dk+ 



(8) 



(2tt) 3 J 2k+ 



J dfc T ^( v / 2fc+) 1 / 2 



x < e 



i (k x —Ict-xt) 



w(s) b(k + , fc-r; s; x + ) 



The factor ^7~7 + here serves to project out the 
components of the quark field ipo that are the in- 
dependent dynamical operators in null plane field 
theory. 

2.3. The quark distribution function 

We are now ready to define the (unrenormal- 
ized) quark distribution function. Let \P) be the 
state vector for a hadron of type A carrying mo- 
mentum P M . Take the hadron to be spinless in 



order to simplify the notation. Construct the 
unrenormalized distribution function for finding 
quarks of flavor j in hadron A as 

f^ A (x)x(P+ / ,P! r \P+,P T ) = 



2x(2tt) 



dk T J2(P + ',PT\ 



xbj(xP + ,kr; s; x ) b 3 (xP + , kr; s; x + ) 



:\P+,P T 



(9) 



In Eq. there are factors relating to the normal- 
ization of the states and the creation/destruction 
operators. There is a quark number operator b^b 
for flavor j. We integrate over the quark trans- 
verse momentum hr and sum over the quark spin 
s. 

2.4. Translation to coordinate space 

With a little algebra, we find 

= ^ / dy-e-™ p+ y- (10) 

x (P+ T |? j (0,1/-, T )7 + V>o j (0, 0, T ) |P+, T ) . 

Notice that the (still unrenormalized) quark dis- 
tribution function is an expectation value in the 
hadron state of a certain operator. The opera- 
tor is not local but "bifocal." The two points, 
(0, y~, Or) and (0, 0, 0t), at which the field oper- 
ators are evaluated are light-like separated. The 
formula directs us to integrate over y~ with the 
right factor so that we annihilate a quark with 
plus momentum xP + . 

2.5. Gauge invariance 

Before turning to the renormalization of the op- 
erator that occurs in the quark distribution func- 
tion, we note that the definition as it stands relies 
on the gluon potential A^(x) being in the gauge 
A + = 0. Let us modify the formula so that 

• the operator is gauge invariant and 

• we match the previous definition in 
A + = gauge. 

The gauge invariant definition is 

/JW^ Jdy~e-™ p+ y- (P+0 T | 



x|P+0 T ) , 



(11) 
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where 

O = Vcxp (^9o£ dz~ A+ a (0,z-,0 T )ta \ .(12) 

Here V denotes a path-ordered product, while the 
t a are the generators for the 3 representation of 
SU(3). There is an implied sum over the color 
index a. 

2.6. Interpretation of the eikonal gauge op- 
erator 

The appearance of the operator O, Eq. (O), in 
the definition ( pi] ) seems to be just a technicality. 
However, this operator has a physical interpreta- 
tion that is of some importance. Let us write this 
operator in the form 

O = Vex.p(-ig J dz~ A+JO, z~, T ) t a j 

xVcxp (itoj dz~ A+ a (0,z- 7 6 T )t a ^ .(13) 

Inserting this form in the definition (|ll]), we can 
introduce a sum over states |iV)(7V| between the 
two exponentials in Eq. (|l^) . We take these states 
to represent the final states after the quark has 
been "measured." 

Consider now a deeply inelastic scattering ex- 
periment that is used to determine the quark dis- 
tribution. The experiment doesn't just annihi- 
late the quark's color. In a suitable coordinate 
system, a quark moving in the plus direction is 
struck and exits to infinity with almost the speed 
of light in the minus direction, as illustrated in 
Fig. ||. As it goes, the struck quark interacts with 
the gluon field of the hadron. 

We can now see that the role of the operator O 
is to replace the struck quark with a fixed color 
charge that moves along a light-like line in the 
minus-direction, mimicking the motion of the ac- 
tual struck quark in a real experiment. 

2.7. Renormalization 

We now discuss the renormalization of the op- 
erator products in the definition (pd|). We use 
MS renormalized fields ip(x) and A^(x) and we 
use the MS renormalized coupling g. The field 
operators are evaluated at points separated by 




Figure 6. Effect of the eikonal gauge operator. 



will be ultraviolet divergences from the operator 
products. We elect to renormalize the operator 
products with the MS scheme. 

For instance, Fig. [?] illustrates one of the di- 
agrams for the distribution of quarks in a pro- 
ton. Before it is measured, the quark emits a 
gluon into the final state. There is a loop in- 
tegration over the minus and transverse compo- 
nents of the measured quark's momentum. This 
loop integration is ultraviolet divergent. To apply 
MS renormalization, we perform the integration 
in 4— 2e dimensions, including a factor (fi 2 e 7 / 4ir) £ 
that keeps the dimension constant while supply- 
ing some conventional factors. The integral will 
consist of a pole term proportional to 1/e plus 
terms that are finite as e — » 0. We simply sub- 
tract the pole term. Notice that MS renormaliza- 
tion introduces a scale /i. 

The definition of the renormalized quark dis- 
tribution function is thus 



dy-er lxP+l 



{p + M 



Ax with Ax^Ax b 



0. For this reason, there 



x^(0,y~,0 T )7 + O ^(0,0, Or) 

x I P+ >0t>ms, ( 14 ) 
where the MS denotes the renormalization pre- 
scription and where 

O =7'exp UoJ^ dz~ A+(0,z-,d T )t a j . (15) 



G 




Figure 7. Renormalization of an ultraviolet di- 
vergent loop integration. 



2.8. Antiquaries and gluons 

We now have a definition of parton distribu- 
tion functions for quarks. For antiquarks, we use 
charge conjugation to define 
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fj /A (x lf j,) = — / dy 



{P + ,0t\ 



x Tr{ 7 +^ (0, IT, T )O ^ (0, 0, T ) } 



where 



O =Pcxp \-ig ^ dz-A+(0,z-,0 T )tij . (17) 

For gluons we begin with the number operator 
in A + = gauge. Proceeding analogously to the 
quark case, we obtain an expression involving the 
field strength tensor F£ u with color index a: 



1 



dy-e~ ixP+ y (P+0 T \ 



xF a (0,y-,0 T ) +u O ab F b (0,0,0 T ) + 



x|^0 t )ms . 



(18) 



3. Renormalization group 

A change in the scale \l induces a change in 
the parton distribution functions f a / A (x, fj,). The 
change comes from the change in the amount of 
ultraviolet divergence that renormalization is re- 
moving. Since the operators are non-local in y~ , 
the ultraviolet counterterms are integral opera- 
tors in k + or equivalently in momentum fraction 
x. Since the ultraviolet divergences mix quarks 
and gluons, so do the counterterms. 

One finds 



dfi 



f tE p */bWS,<*.0*)) h/A{t,n). (20) 



The Altarelli-Parisi (= GLAP = DGLAP) kernel 
P a /b is expanded in powers of a s . The a], and a 2 s 
terms are known and used. 

3.1. Renormalization group interpretation 

The derivation of the renormalization group 
equation ( p0| ) is rather technical. One should not 
lose sight of its intuitive meaning. Parton split- 
ting is always going on as illustrated in Fig. ||. A 
probe with low resolving power doesn't see this 
splitting. The renormalization parameter \x cor- 
responds to the physical resolving power of the 
probe. At higher /i, field operators representing 
an idealized experiment can resolve the mother 
parton into its daughters. 




ITTTTTT 



where 



= Texp\igJ dz~ A+ (0, z~, T ) t c I . (19) 



Figure 8. A quark can fluctuate into a quark plus 
a gluon in a small space-time volume. 



Here the t c generate the 8 representation of 
SU(3). 
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3.2. Renormalization group result 

One can use the renormalization group equa- 
tion ( po|) to find the parton distributions at a 
scale fi if they are known at a lower scale (1q. 
Fig. ^| shows an example, the gluon distribution 
at (i = 10 GeV and at (X = 100 GeV (using the 
CTEQ3M parton distribution set). Notice that 
with greater resolution, a gluon typically carries 
a smaller momentum fraction x because of split- 
ting. 



0.1 



x 

Cm 

CI 

X 




100 GeV 



0.001 0.01 
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0.1 



Figure 9. Evolution of the gluon distribution be- 
tween fi = 10 GeV and (i = 100 GeV. 



4. Translation to local operators 

We have defined the parton distributions as 
hadron matrix elements of certain operator prod- 
ucts, where the operators are evaluated along a 
light-like line. Now we relate the parton dis- 
tributions to products of operators all at the 
same point. It is these local operator products 
that were originally used in the interpretation 
of deeply inelastic scattering experiments. For 
lattice QCD, evaluation of operator products at 
light-like separations would be, at best, very diffi- 
cult, so the translation to local operator products 
seems essential. 



4.1. Quarks 

Note that, according to the definitions ( |l4| ) and 
©, 

fj/A{x,V) = 0, X>1, 
fj /A (x,fx) = 0, x>l, 
fj/A(~x,n) = -f] /A (x,(i). (21) 

Consider the moments of the quark/ antiquark 
distributions defined by 



71/ 



, (/*) = 
1 dx 



(22) 



for J = 1,2, 
is 



x J {fj/A{x;v) + (-i) J f- j/A (x;iJ.)} 



My (ji) = / —x J f 3/A {x-(i) 



Given the properties above, this 
dx 



(23) 



Obtaining dx is the essential step. 

From the operator definitions of the /s, this is 

M ( /\») = 

1 



dx e - lxP+v ~ 



-i d 



j-i 



v P+ dy- 

*(PU> 3 {y-)i + o{ y -^) i>M\P)m ■ ( 24 ) 

Performing the x-integration gives a 5(y~). Thus 
we get a local operator. The d/dy~ differentiates 
either the quark field or the exponential of gluon 
fields in O, Eq. (jl|). We find 



Mf'(n) = 



1 
2 

where 
D* 1 



J-i 



^•(0)|P) 



dy^ 



(25) 



MS ' 



(26) 



We have now related the moments of the quark 
distribution to products of operators evaluated at 
the same point. However, this is not yet ready 
for the lattice because it would not be easy to 
differentiate the operators with respect to y~ . To 
obtain a more useful expression, consider ) 
defined by 

{pwpft...pw} K = 



(27) 



(P\1>i (0) { 7 M1 ■ ■ ■ iD» } TS ^ (P)\P)m 
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where TS denotes taking the traceless symmetric 
part of the tensor enclosed. 
Then 

(0«)=Mj J V)= (28) 
f 1 dx 

/ — x J {f j/A (x;fi) + {-l) J f- j/A {x;n)} . 
Jo x 

This is our final result. 

We can now imagine the following program for 
a lattice calculation. One could measure (O^) 
on the lattice for J = 1,2, .... Of course, this 
would not be so easy, but it would gives moments 
of the quark/antiquark distributions that could 
be compared to the corresponding moments of 
the distributions determined from experiments. 

4.2. Gluons 

We follow a similar analysis to relate the gluon 
distribution function to local operator products. 
From the definition (|l8|), it follows that 

fg/A{x,H) = 0, X>1, 
fg/A(-X,V) = ~fg/A{x,n). (29) 

For J — 2,4,6,..., we relate moment integrals 
over the interval < x < 1 to moment integrals 
over all x 

f —x J f g /A{x]n) = \ [ —x J f g/A (x;n).(30) 

JO x L J -oo x 

Define operator matrix elements (Og) for J = 
2,4,6, ...by 

{pwpMi . . . pmj} ts p(J)} _ ( 31 ) 

1{P\{F^(0)W^ ■ --iD^ F^(0)} TS \P) m . 

Then these operator matrix elements are the mo- 
ments of the gluon distribution: 

(0«)=Mf>(/i) = f -x J f g/A {x^). (32) 
Jo x 

Again, one can contemplate measuring {O g J ^) on 
the lattice for J = 2,4,.... The result could be 
compared to the corresponding moments of the 
phcnomenological gluon distribution. 



5. Determination from experiments 

I have alluded several times to the fact that 
parton distribution functions are related to ex- 
perimental results. Let us look at this in some 
more detail. 

5.1. Example 

Fig. |l0| shows data for the cross section to ob- 
serve one jet plus anything in high energy pp col- 
lisions. The data from the CDF Collaboration 
H are compared to theory using CTEQ3M par- 
ton distributions that did not have this data as 
input. (In the figure, Et denotes essentially the 
transverse momentum of the jet.) 
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Figure 10. One jet inclusive jet cross section 
da/dEx, (Data — Theory) /Theory versus Et- 
Source: CDF. 



One notes that the data do not agree with the 
theory at the highest values of Et- Does this 
indicate a breakdown of QCD at the smallest dis- 
tances? An alternative explanation is that 
the gluon distribution used in the theory was ac- 
tually poorly constrained at large x (x ~ 0.5) 
by other experiments. In support of this expla- 
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nation, the CTEQ collaboration offered the the- 
ory curve shown in Fig. [llj This theory curve 
was obtained using modified parton distributions 
that still fit other data well. (Note: the data in 
this figure predates that in Fig. |Io] and has lower 
statistics.) 



Table 1 

Moments of gluon distribution for two sample 
parton distributions sets 

CTEQ3M CTEQ4HJ 

Mf ] (2 GeV) 4.5 x 1CT 1 4.3 x 1CT 1 

Mf ) {2 GeV) 1.6 x 1CT 2 1.9 x 1CT 2 

M ( g 6) (2 GeV) 2.1 x 1CT 3 4.4 x 1CT 3 



CDF Inclusive Jet Cross Section (Preliminary) 
NLO OCD using CTEQ3M, ^=Pt/2 



Norm= 1.0 Jet-Fit 

Latest CTEQ Fit (Norm = 0.93) 



Correlated Systematic Uncertainties 

-0.4 - 

" i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i ~ 
0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 



Figure 11. (Data — Theory)/Theory versus xt = 
2E T /V* for CTEQ4HJ partons, "Norm = 1.0 
Jet-Fit." Source: CTEQ. 



Could a lattice measurement of moments of the 
gluon distribution distinguish between these pos- 
sibilities? In Table [l], I show the lowest moments 
of the gluon distributions from the two parton 
distribution sets in question. The two gluon dis- 
tributions are very different at large x, but almost 
the same at small x. Thus one has to go to the 
J = 6 moment to see a substantial difference. 
Measurement of M^ 6) (2 GeV) in a lattice calcu- 
lation would, I think, be quite difficult. 

5.2. What parton people do 

There are three main groups currently doing 
"global fitting" of parton distributions: 

• Martin, Roberts and Stirling (MRS) (n), 

• Gluck, Reya, and Vogt (GRV) pf, 



• Botts, Huston, Lai, Morfin, Owens, Qiu, 
Tung, Weerts,. . . (CTEQ) §. 
Other groups have done this in the past, e.g. 
Duke and Owens and Eichten, Hinchliffe, Lane 
and Quigg. 

The idea of the global fitting is to adjust the 
parton distribution functions to make theory and 
experiment agree for a wide range of processes. 
For example, recent CTEQ fits have used the fol- 
lowing processes: 

-p -> X 

- Fe -> X v 
p + p->fi + ju + X p + 2 H ->■ /i + p + X 
p + Cu — ► n + p, + X 

p + p^ w ^i + i+x 

p + p — > 7 + X . 

The main features of a program of global fit- 
ting are as follows. One chooses a starting 
scale no (say 2 GeV ). Then one writes the 
fa/ P {x]Ho) m terms of several parameters for 
a = g,u, u, d, d, s, s. (Typically the heavy quark 
distributions, for a = c, c, b, b, t, t, are generated 
from evolution, not fit to data.) For example, one 
may choose 



e - 

M ' 
v ■ 



- 2 H 
-Fe 
2 H- 



> X 

> X 



f(x; f i () )^Ax B (l-xf(l + Ex D ) 



(33) 



The parton distribution functions obey certain 
flavor and momentum sum rules, such as 



dx [f u / p (x,n) -f u /p(x,n)] = 2, 

Y] / dx xf a/p (x,fi) = 1- (34) 
„ Jo 



These sum rules follow from the operator def- 
initions (0), (Jig), ((HI). The parametrizations 
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are chosen so that the flavor and momentum sum 
rules are obeyed exactly. 

Now one picks a trial set of parameters. This 
determines f(x;/j,o), from which one calculates 
f(x;/j,) for all \i by evolution. Next, given the 
f(x; n), one generates theory curves for each type 
of experiment used. Finally, one compares the re- 
sults to the data. 

This sequence is iterated, adjusting the param- 
eters to get a good fit. There are more than 1000 
data and only about 25 parameters to be fit, so 
the fact that this procedure works at all is an in- 
dication that QCD is the correct theory of the 
strong interactions. 

6. Relation to lattice QCD 

Let me give here a personal view about how 
lattice calculations of the lowest moments of par- 
ton distributions, continuing the work reported 
at this conference |T^ |, might enhance our under- 
standing of particle physics. 

First, what do we know now? To a reasonable 
approximation, the parton distribution functions 
have already been determined from experiment. 
Furthermore, there is partial understanding of 
parton distributions in two limits: one has some 
idea about x — > behavior from Regge physics 
and from the evolution equation; one has some 
idea about x — > 1 behavior from Feynman graphs. 
However, this understanding is of a rather limited 
nature. In general there is little theoretical un- 
derstanding of why parton distributions are what 
they are. 

A lattice calculation of moments of parton dis- 
tributions could, first of all, test QCD. One could 
calculate something about hadron structure in- 
stead of simply measuring that structure. It is 
possible, in principle, that lattice calculations 
could add to the precision of parton determi- 
nations, particularly for the gluon distribution. 
However it seems to me that the goal of doing 
better by calculation than by experiment is rather 
distant at this point. Finally, and perhaps most 
importantly, it seems to me that lattice calcula- 
tions might lead to an improved understanding 
of hadron structure by either validating or falsi- 
fying simple pictures about how the partons are 



arranged in a hadron. I turn to this question in 
the sections that follow. 

7. Are partons inside "valence quarks"? 

Figure [l2| illustrates a simple picture of how 
partons might be arranged in a proton. A proton 
consists of three valence quarks, U U D. A valence 
quark in this sense is a big object. Each valence 
quark consists of short distance quarks and glu- 
ons. I do not know the history of this simple 
idea, but one illustrative paper on the subject is 
Ref. @. 




Figure 12. A hadron may contain valence quarks 
that in turn are made of short distance quarks 
and gluons. 



This idea may be expressed in the formula 

fa/A{x\n)= [ f £ F J/A (Ofa/j(x/C,ri-(35) 
Jx *> J=U,D 

Here Fj/a(£) represents the distribution of va- 
lence quarks of flavor J in a hadron of type A. 
Then f a /j(x/^; fx) d(x/() is the probability to find 
a parton of type a carrying a fraction x/£ of the 
valence quark's momentum. 

Note that, without an additional hypothesis 
for Fj/a(£), this formula has negative predic- 
tive value: in order to predict the values of 
the functions f a / P {x;n) we introduce unknown 
functions f a /j( x /^^) an d Fj/a(Q- However, 
if you look not just at hadrons p and n but 
at hadrons p, n, n, p,u>, . . . then there is predic- 
tive power. There is experimental information 
on the distribution of hadrons in the pion from 
7r + p — > fi + + fi~ + X. In addition, there is 



some experimental information on the distribu- 
tion of hadrons in the p° from deeply inelastic 
scattering from photons. However, experimen- 
tal information on the distribution of partons in 
other hadrons is lacking. Here one sees an advan- 
tage of calculation: it is much easier to create all 
kinds of hadrons on a lattice than in an accelera- 
tor. 

7.1. Possibilities for Fj/a(x) 

One simple hypothesis for Fj/a(x) would be 
that a constituent quark is simply an elemen- 
tary quark as seen at a very low resolution 
scale (i. Thus one might identify Fj/^(x) = 
fj/A ( x > 0.3 GeV). However, this hypothesis does 
not seem to work jl^] . 

Another possibility is that the constituent 
quarks are very weakly bound in a hadron. Then 

F u/p {x) » 2 5(3!- 1/3), 

F D/p (x) « 15(3! -1/3). (36) 

If one investigates only protons and neutrons, 
then this hypothesis has zero predictive power. 
(That is better, of course, than negative predic- 
tive power.) However, of one looks at p, n, p, to... 
then there is some predictive power. 

7.2. Parton correlations 

The valence quark picture appears to say some- 
thing about multiparton correlations. Examining 
Fig. flll|), we see that the partons in a proton 
are clumped together more than they would be if 
they were spread throughout the proton. Thus a 
suitably defined two parton correlation function 
would be larger at short distances than it would 
be in an uncorrelated model. 

There is some experimental information about 
two parton correlations from double parton scat- 
tering. Such an experiment is illustrated in 
Fig. || 

One could investigate this issue in a lattice cal- 
culation. One would define a two parton distri- 
bution function f( a ,b)/A( x a, %b, r) where x a and Xb 
are the momentum fractions of two partons and 
r is the transverse separation between them. The 
simplest uncorrelated model (ignoring even cor- 
relations due to momentum conservation) would 
be that this function is the product of f a /Ai x a)i 
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Figure 13. If the short distance partons are con- 
tained inside valence quarks, then the probability 
that two partons are close together is larger. 



fb/A{ x b), and a theta function that requires that 
both partons are not more than a distance R from 
each other, 9(r < R). There must be more corre- 
lation than this. The question is, how much? 

8. Partons in the pion cloud 

Figure [l^ illustrates another simple picture, 
usually known as the Sullivan model [T^ ]. Some- 
times a proton is a proton. Sometimes a proton 
is a proton plus a pion. More generally, a pro- 
ton can exist virtually as a baryon plus a meson. 
Each baryon or meson consists of short distance 
quarks and gluons. 

A formula for the parton distributions in this 
model is 

fo/A(x;n)= f ^-J2 F h/a(0 L/h(x/^^.(37) 

J x S 

Here H represents hadrons that might be virtual 
constituents of the proton and F H / A (£) represents 
the distribution of these virtual hadrons in the 
proton. One makes a model for Fh/a{0 based 
on measured meson-baryon couplings. 

Note that A could be a proton, a "physical pro- 
ton," while H could be a proton, a "bare proton." 
The relation between these is, in my opinion at 
least, rather murky. Nevertheless, there must be 
some truth to the model. If there is some truth to 
it, the model suggests that sometimes a hadron 
fluctuates to a much bigger size. Again, an inves- 
tigation of multiparton correlations could help to 
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Figure 14. Parton correlations can be examined 
experimentally by looking for double parton scat- 
tering. 





Figure 15. A physical hadron may consist of bare 
hadrons that in turn are made of quarks and glu- 
ons. 



pin the question down. 

9. Conclusion 

I conclude that lattice calculations continuing 
the work [|l3| presented at this conference could 
test QCD by predicting parton distributions di- 
rectly from the QCD lagrangian. One also has 
the possibility of asking and answering questions 
about hadron structure that go beyond what can 
be easily measured experimentally. 
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